Universal energy distribution for interfaces in a random field environment 
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We study the energy distribution function p(E) for interfaces in a random field environment at 
zero temperature by summing the leading terms in the perturbation expansion of p(E) in powers 
of the disorder strength, and by taking into account the non perturbational effects of the disorder 
using the functional renormalization group. We have found that the average and the variance of 
the energy for one-dimensional interface of length L behave as, (E)r oc LlnL, AEr oc L, while 
the distribution function of the energy tends for large L to the Gumbel distribution of the extreme 
value statistics. 

PACS numbers: 05.20.-y, 74.25. Qt, 64.60.Ak, 75.60.Ch 



I 

C/3 



I 

C 

o 
o 



CN 
> 

-sj- 
CN 
O 

CO 

O 

03 



i 

a 
o 

o 



X 



The concept of energy landscapes is of current inter- 
est in different topics such as structural glasses, spin 
glasses, proteins, flux lines etc. [213 ■ The existence 
of metastable states is crucial for the complex behav- 
ior in these systems. The domain wall counterpart of the 
random field Ising model provides an example of the 
problem which possesses complex properties, and can be 
quantitatively treated using the well established analyti- 
cal methods such as the functional renormalization group 
(FRG) method [T(j and the method of replica symme- 
try breaking (RSB) [13. A significant progress has been 
achieved in recent years in understanding the behavior 
of interfaces in disordered media at equilibrium llOl and 
the driven interfaces at the depinning transition [12h14| . 
It is expected that in equilibrium or below the depinning 
transition there are many metastable states constituting 
the energy landscape. This makes the interface problem 
a natural candidate to study the concepts of energy land- 
scapes. For recent theoretical and numerical studies of 
the related systems under the perspective of the energy 
landscape see Refs. [l5ll2l| . In this Letter we present the 
results of the study of the distribution function (DF) of 
the energy p{E) for an interface in a random field envi- 
ronment at equilibrium. The main result of this Letter 
is that for large interfaces p(E) is a universal function 
which coincides with the Gumbel distribution of the ex- 
treme value statistics. The dynamic formalism we use 
here can be applied to the study of the dynamic aspects 
of energy landscapes such as the relaxation of the energy, 
the two times energy correlation functions etc. 

The interface motion in a disordered medium at T = 
is described by the equation 



_! dz(x, t) 
dt 



= 7V 2 z + F + g(x,z), 



(1) 



where fi is the mobility, 7 is the stiffness constant, and F 
is the driving force density. The quenched random force 
g(x, z) is assumed to be Gaussian distributed with the 
zero mean and the second cumulant (g(x, z)g(x', z')) = 
5 d (x~x')A(z — z'), where A(z) is a short ranged function 
with the width a, and d is the interface dimension. 



It is well-known that the Langevin equation JQ) can be 
reformulated in terms of the Fokker-Planck equation for 
the conditional probability density P(z(x), t; z°(x), t°) to 
have the profile z{x) at time t by having the profile z (x) 
at time t°. This Fokker-Planck equation can be written 
as an integral equation, which, for an interface of a finite 
length L, reads 

P(z,t;z°,t°) = P (z,t;z°,t ) 

- fi J dt' J T>z'P (z,t;z',t') 

x Y, d *' k ,9k>{z'W,t'\z°,t Q ), (2) 



where Zk = J d d xz(x) exp(-ikx) and gk{z) — 
J d d x exp(— ikx)g(x, z) (k = (fci, kd), h — 2nji/L, 
ji = 0, ±1, ...) are the Fourier transforms of the inter- 
face height and the quenched force, respectively. J T>z in 
^ stays for integrations over the modes z = {zfc}. The 
bare conditional probability for non zero modes reads 

P (z,t;z°,t a ) = l[S(z k - z° k exp(- lf ik 2 (t - t ))). (3) 

k 

Analogously to the case of one Brownian particle the 
formal solution of Eq. can be represented as a path 
integral. 

The probability DF of the energy E(z) = E c \{z) + 

Edis{z) = / d d x(^(\7z) 2 — g(x, z)) can be calculated 
using the conditional probability density P(z(x), t; 0, 0) 
as follows 

p(E(t)) = J Vz(x)S(E - E(z))P(z(x),t;0,0). (4) 

It is convenient instead of p(E(t)) to consider its Fourier 
transform p(s) which is obtained from Eq. as 



n—0 ' rn—0 



(5) 



In this Letter we will restrict ourselves to the study of 
the energy DF in the steady state, i.e. for t — > 00. In 



2 



this limit {E™E^ is m ) is related to the static equilibrium 
correlation function 

(z(x\)z{x2)) = lhn / T>z(x)z(x\)z(x2) 
*— 00 J 

xP(z(x) t t;0,0) = J T^^e*^-^). (6) 

Let us for simplicity elucidate the computation of the 
nth moment of the elastic energy 



<33> = (~) T 



Vz / k\\z kl \ 



kl\z kn \ 2 POM; 0,0), 
(7) 

where we expressed _E c i through the Fourier components 
of the interface height z{x). For an interface of a finite 
size L the integral f. means L~ d J2 k - To compute (J7J to 
the lowest order in disorder strength we iterate Eq. (J2J 
2n times and insert it into (JJJ. Expecting that the steady 
state does not depend on the initial interface configura- 
tion we have taken the latter in © and J7J) to be flat 
at to — 0. The average over the random forces, which is 
carried out by using the Wick theorem, yields connected 
and disconnected expressions. The connected expression 
contains only one integration over k, while the number of 
integrations over k in a disconnected expression is equal 
to the number of connected parts in that expression. Let 
us consider the calculation of the connected part of (E^). 
As a result of integrations by parts in with P(z, t; 0, 0) 
being iterated 2n times the 2n derivatives with respect 
to zL (sec Eq. |J5J) will act on z ki in J7J. This has the 

consequence that pairs of 2n momenta k[, k' 2n asso- 
ciated with the right-hand side of (being iterated) 
become equal consecutively to one of k\, k n in J7J. 
The number of such possibilities is (2n)\. The rid of 2n 
ordered time integrations in P{z, t; 0, 0) gives the factor 
l/(2n)! The number of possibilities to get a connected 
loop diagram shown in Fig. ^ with n continuous lines is 
2"~ 1 (n— 1)!. Integrations over xx, X2n-i arising from 
the above expression of g k (z) provides that the momenta 
of the modes being connected by a dashed line, which 
is associated with the disorder correlator, become equal. 
The integration over xi n gives the factor L d . The inter- 
mediate z' k are zero for flat initial interface configuration 
due to delta functions in Q. As a result the arguments 
of disorder correlators A(z) become zero. Collecting all 
combinatorial factors and taking the limit t — > 00 we find 
the following expression of the connected part of (E^) 



>a = ^A ( or 7 -^/ fe ^. 



(8) 



The computation of (E^E2 ia m ) to the same order 
is similar and gives {E£E™r m ) = i(-2)"" m (n - 
iy.A(0) n -f- n L d j k l/(fc 2 )™. The use of the latter yields 
(E n ) c which is obtained from (JHJ multiplied by the 
factor (—1)™. The expression A (E n ) c is associated with 
the loop diagram consisting of n continuous lines (see 



+ ... 



FIG. 1: The loop expansion of connected diagrams contribut- 
ing to the energy distribution function. 



Fig. [JJ . The factor 2n is the symmetry number of the 
diagram. The straightforward analysis gives that the ex- 
pansion JSJl can be represented as a series of loop dia- 
grams. The use of the connectivity theorem enables us 
to write the Fourier transform of the energy DF (JjjJ as ex- 
ponential of the series of connected loop diagrams shown 
in Fig. [2 



Pis) 



exp 



is A(0)\ 



(9) 



Notice that p(s) given by the diagram series in Fig. 
is closely related to the loop expansion of the effec- 
tive potential in Quantum Field Theory studied in (23j. 
Replacing the integral in by the sum according to 
L f k f(k) - ET=-oo f(2nj/L) we find in d = 1 



p(s) = J] (1 + isE /f 
i=i 



7TV isE 
smh(n^/IsE~o) ' 



(10) 



where Eq = — A(0) L 2 / (An 2 ^) is the characteristic energy 
for an interface with the perturbational roughness w cx 
i 3 / 2 , which follows from w cx L^ d ^ 2 for d = 1. 

Eq. (|10f) has only simple poles s — ij 2 /Eo in the lower 
half-plane, so that the inverse Fourier transformation of 
(|10fl can be easily performed as a sum over all poles by 
using the Jordan's lemma. As a result we obtain the DF 
as p{E) = \E \- X f (E/E ), E < , where 



i=i 



1) 



U- 



(11) 



The derivation of p(E) shows that 111 1[1 describes also 
the DF of the elastic energy [x — E/\E \ > 0) and 
the disorder energy {x — E/2E > 0). Eq. lfTT|) co- 
incides with the dimensionless width DF for the one- 
dimensional random- walk interface studied in Ref. |24| . 
Using the method of stationary phase it was shown in 
24] that JTTJ can be well approximated for small x by 
f(x) w v ^Jx~5(n 2 /2 - x)e^ 2 / ix . Using JTTJ we have 
computed the average energy, (E) = n 2 E /6, and the 
variance AE = {(E 2 ) - {E) 2 ) 1 / 2 = tt 2 | 1 / (3^X0) . 
Eq. i|ll|) is the exact perturbational result generaliz- 
ing the result established by Efetov and Larkin [2{| for 
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the height-height correlation functions JHJ) , which can be 
readily proved using supersymmetry [2q |. Nevertheless, 
both © and lfl"T|l are wrong due to the fact that © gives 
the value (4 — d)/2 for the roughness exponent instead of 
the correct Imry-Ma value ( = (4 - d)/3. 

We now will take into account the effect of the renor- 
malization on the e nerg y DF using the results of the FRG 
|lOj| . According to |lOj the renormalized disorder corre- 
lator Ar(0) acquires the scale dependence i 2 '- -6 . Taking 
into account the latter in © gives the correct value of 
the roughness exponent To enable a crossover to the 
perturbational regime at small length scales we use the 
ansatz 



A fl (0)=A(0) l + (k/k c ) 



\*C- 



(12) 



where e = 4 — d, and the wave- vector k c = 2tt/L c is asso- 
ciated with the Larkin length L c ~ (-/ 2 a 2 / A(0)) 1/e . The 
ansatz ()12|l describes the scale dependence of A/j(0) at 
the cusped fixed-point solution of the disorder correlator, 
A R (0) ~ A(0) [k/k c } £ ~ 2C for k < k c , and describes the 
crossover to the perturbational regime, Ar(0) ~ A(0) 
for k ^S> k c . Using the renormalized Afl(O) in Eq. © we 
obtain the Fourier transform of the renormalized distri- 
bution of the energy in d = 1 as 



Pr(s) 




isEo 



(13) 



where r/ = L/L c and E Q = — A(0) L C L/ (in 2 j). Carrying 
out the inverse Fourier transformation of (|13J) we obtain 
p R {E) = \E \^ 1 f R (E/E ; rj), where the function f R (x; 77) 
is given by 

f R (x-n) - V( 1V+1 ± 11 ± + 2 3h) c -m+iMx 
TR{x,V)-^ L > r(n + l)r(7) 



3=1 



(14) 

For string lengths much shorter than the Larkin length, 
77 <C 1, the DF (|14f> pass over to the perturbational result 
(|ll|l . Similar to the height-height correlation function at 
equilibrium we expect that (|14|) . which is the result of 
the renormalization of (|11(1 to order e, is exact. Eq. Ijl4|l 
applies to order e at the depinning transition too with 
the difference that in this case there are corrections to 
(|14fl of order e 2 . However, we expect that the latter will 
be small as it is the case for corrections of order e 2 to the 
interface width distribution at the depinning transition 
[2l| . The average energy (E) R derived from Eq. (|T3^I is 



[V(r) + 1) + C] E 



~ [lnri + C]E + O(l/ri) (x LI11L, (15) 

where ^(x) is the digamma function and C — 0.5772... 
is the Euler's constant. The energy fluctuation AE R ob- 



tained from (|13fl reads 



AE R = I S 1 



— + &(r l + l)-2(C + $(ri + l))/j ] 



nl/2 



— |£ |+O (In V / V ) cxi. 



(16) 



The result, AE R oc L, agrees with the estimate of the 
energy by using dimensionality arguments with correct 
roughness exponent £. Notice that due to the logarithmic 
term in (|15|1 (E) R and AEr scale in different way, and 
the relative fluctuation, AEr/ (E) r, disappears as 1/ InL 
for large L, which is in contrast to 1/VZ behavior for a 
Gaussian distribution. The latter reflects the relevance of 
fluctuations over all length scales. The higher moments of 
the energy distribution {HJl scale as ((E- (E)) n ) oc AE R . 
Fig. |2 shows fR(x;rj), which is given by Eq. I|14() . as a 
function of x — In rj. 




FIG. 2: The renormalized distribution of the energy for a line 



in a random field environment. Dashed line: L/L c 
solid line: the Gumbel distribution. 
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We now will consider the asymptotic behavior of i|14|) 
in the limit of long lines, L L c . Changing x in favor 
of x — In rj = y and taking the limit 77 — > 00 we calculate 
the sum over j in (|14fl and arrive at 



fn{y) = ~P(y) = exp(-y - exp(-y)), 



(17) 



which is nothing but the Gumbel distribution of the ex- 
treme value statistics |2^. The universality of /i?(y) is 
due to the universal character of fluctuations on large 
scales, which are described by the fixed-point solution of 
the FRG 01 • Notice that the expectation value of y cal- 
culated with l|17|l gives the Euler's constant C which is 
in consistence with Eq. (|15fl of the average energy. We 
have checked that the limit of the distribution Jr{x; rj) 
for 77 — > 00 is insensitive to the details of the renormal- 
ization at scales smaller than the Larkin scale. 

The Gumbel distribution is one of the three possible 
limit distributions in the extreme value statistics |29| . 
which concern the distribution of the maximum M n = 
max{£i, £„} of the set of identically distributed ran- 
dom variables £j (i = 1,2, ...,n). The asymptotic distri- 
bution P n (x) for M n in limit n — > 00 does not depend on 
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details of the distribution of £j and under fulfilling some 
conditions [2!j has the form P n (x) ~ V(x — Inn) where 
■P(y) is given by Eq. (|1T|) . The combination y = x — Inn, 
where n is the number of random variables guarantees 
that the distribution remains invariant for n — > oo. 

Vinokur ef a/. |lfij ] have used the Gumbel distribu- 
tion to describe in a phenomenological way the energy 
barriers distribution for a flux line in a random environ- 
ment. The creep motion of the flux line in the limit of 
small driving force F and low temperature is controlled 
by thermally activated jumps. The thermally activated 
advance of the flux segment of length L is controlled by 
the global barrier U = max{t/i, ...,[/„}, where Ui is the 
barrier for the subsegment i of length L c with the num- 
ber of subsegments n = L/L c . It was suggested in 0] 
that the probability distribution for a given segment L is 
V(U/U C - \n(L/L c )), where U c ~ ja 2 L*- 2 is the mini- 
mum average barrier between neighboring metastable po- 
sitions of a pinned segment L c , so that the typical barrier 
of a segment of length L scales then as U oc U c ln(L/L c ). 
Bouchaud and Mezard [l^ showed that the Gumbel dis- 
tribution describes the energy distribution in a class of 
random energy models possessing the one-step RSB. The 
Gumbel and related distributions were used in [3(|, to 
describe universal fluctuations in correlated systems. It 
was shown in |3l| that the Gumbel distribution appears 
in systems with 1// power specta. 

The relation of the energy distribution with the ex- 
treme value statistics can be illustrated as follows. It is 



known from the treatment of the problem in the frame- 
work of the replica variational approach |llj| that the 
system under consideration demonstrates RSB |22| ■ The 
RSB is related to the existence of many metastable states 
and to the lost of the ergodicity of the system. On the 
contrary, the perturbational result l|llfl corresponds to 
the average over states which belong to different noner- 
godic subsystems. The metastable states, which corre- 
spond to the minima of the total energy of the interface, 
are rare states among other states of the energy land- 
scape, and has to reflect the statistics of these extreme 
states . It was shown in 0, 0] that energy minima 
behave in the same way as barriers, which as suggested 
in are distributed in accordance with the Gumbel 
distribution. 

In conclusion, we have studied the distribution func- 
tion of the elastic, disorder, and the total energies of an 
interface in a random field environment by summing the 
leading terms of the perturbation expansion in powers of 
the disorder. The nonperturbational effects of the dis- 
order are taken into account using the FRG method. 
We have found that the average and the fluctuation 
of the energy for one-dimensional interfaces behave as, 
(E)r oc L\nL, AEr oc L, while the energy DF tends for 
large L to a universal function which coincides with the 
Gumbel distribution. 
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